Crystalline and Electronic Structures of Molecular Solid C50CI10: First-Principles Calculation 
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A molecular solid C50CI10 with possible crystalline structures, including the hexagonal-close-packed (hep) 
phase, the face-centered cubic (fee) phase, and a hexagonal monolayer, is predicted in terms of first-principles 
calculation within the density functional theory. The stable structures are determined from the total-energy 
calculations, where the hep phase is uncovered more stable than the fee phase and the hexagonal monolayer in 
energy per molecule. The energy bands and density of states for hep and fee C50CI10 are presented. The results 
show that C50CI10 molecules can form either a hep or fee indirect-gap band insulator or an insulating hexagonal 
monolayer. 
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PACS numbers: 71.20.Tx, 31.70.Ks, 73.61.Wp 

Fullerenes such as Ceo, C70 and their larger homologs sat- 
isfy the so-called isolated pentagon rule (IPR)Clli. This rule 
states that the most stable fullerenes are those in which every 
pentagon is surrounded by five hexagons. For some time, peo- 
ple believe that the non-IPR fullerenes could be unstable due 
to their adjacent pentagons, and may be difficult to synthesize. 
In recent years, it has been uncovered, however, that the IPR 
can be violated by metallic endohedral or exohedral chemical 
deriving fullerenes||^|4^(J,0,0,|3l. On account of the adjacent 
pentagons and the high curvature of molecular surface, the 
non-IPR fullerenes might have unusual electronic, magnetic, 
and mechanical properties. Last year, Xie et al.\9] had suc- 
cessfully synthesized a new exohedral chemical deriving non- 
IPR Ds/i fullerene[50], say, C50CI10, with 10 chlorine atoms 
added to the pentagon-pentagon vertex fusions. The molecu- 
lar structure of C50CI10 looks like a Saturn-shaped profile, as 
depicted in Fig. 1, which is compatible with the Euler theo- 
rem, ensuring that the fullerenes with fewer than sixty C atoms 
and comprised of only pentagons and hexagons do not satisfy 
the IPR. The existence of D^h C50CI10 is evidenced by mass 
spectrum, and ^^C NMR spectra; and of particular interest 
is that C50CI10 could easily react with a variety of organic 
groups to form new compounds! 9| which might show inter- 
esting chemical and physical properties. Since its discovery, 
there have been several investigations LI (ill 1[I12lI13[Ii4I15Ii 
on the electronic and optical properties of C50CI10 clusters. 

So far, none has discussed, however, the possibility if 
C50CI10 molecules can form a solid. As is well-known, Cgo 
clusters can condense to form a solid with either the face- 
centered-cubic (fee) phase| 16] or the hexagonal-close-packed 
(hep) phaseQTJ. Although the molecular structure is different, 
there is no reason to believe that C50CI10 molecules cannot 
condense to form a solid. Therefore, it would be interesting 
to address this issue at least from a viewpoint of numerical 
simulations. In this paper, we report a first-principles study 
on the possible solid CsoClin within the framework of density 
functional theory (DFT) |18] with local density approxima- 
tion (LDA)1 19|. From the total-energy calculations we have 
found that C50CI10 molecules could form a stable hep or fee 
crystal or a hexagonal monolayer, though the hep phase is 
more stable than the fee phase and the hexagonal monolayer 
in energy per molecule. The optimized lattice constants are 




FIG. 1: (Color online) Schematic structure of molecule C50CI10: (a) 
top view, and (b) side view. The dark grey and green (light grey) 
balls represent C and CI atoms, respectively. The numbers I to IV 
indicate the four types of C atoms. 



ao = 13.57 A, cq = 11.39 A for the hep phase, qq = 19.20 
A for the fee phase, and ao = 13.52 A for the hexagonal 
monolayer The energy bands and density of states (DOS) 
for three structures have been calculated. The results show 
that the solid C50CI10 is a band insulator with an indirect en- 
ergy gap of 1.59 eV for the hep phase or 1.79 eV for the fee 
phase, and with a gap of 1.85 eV for the hexagonal mono- 
layer, where the gap widths are comparable with those of Ceo 

solids m. 

The total energies, energy bands and DOS are calculated 
by means of the ABINIT package 1 22]. This package is 
coded within the DFT framework based on pseudopotentials 
and plane waves, which relies on an efficient fast Fourier 
transform algorithm l!2^ for the conversion of wave func- 
tions between real and reciprocal spaces, on the adaptation 
to a fixed potential of the band-by-band conjugate gradient 
method! 24|, and on a potential-based conjugate-gradient al- 
gorithm for the determination of the self-consistent potential 
(2,^ El- Troullier-Martins norm conserving pseudopo- 
tentials 1 28] are applied to mimic the electron-ion interac- 
tion, and the Teter parametrizationf25] of the Ceperley- Alder 
exchange-correlation potential is used. The kinetic energy 
cutoff in the plane-wave basis is taken as 20 Hartree, and the 
tolerance for absolute differences of the total energy is set as 
10^^ Hartree. The convergence of the total energy to the ki- 
netic energy cutoff has been checked. 

As shown in Fig. 1, the C50CI10 molecule has 4 unique 
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FIG. 2: The presumed crystalline structures: (a) a tetragonal lattice; 
(b) a simple hexagonal lattice; and (c) a hep lattice. Note that to show 
lattices clearly, the C50CI10 molecules are drawn much smaller than 
those in the real proportion. 



types of carbon atoms (I to IV) within D^h symmetry, and 
ten CI atoms are added to the most reactive Cjv sites (i.e. 
to pentagon-pentagon vertex fusions). Such a structure obvi- 
ously violates the IPR. The C-Cl bond length in C50CI10 is 
1.80 Aj^, which is comparable with that in CCI4, 1.78 A. 
We have calculated the valence-electron density of C50CI10, 
and found that the high density of the C-Cl bond shows that 
the ten CI atoms are indeed bonded with the C50 cage, and a 
stable C50CI10 molecule can be formed. 

To seek for stable solid phases, by carefully analyzing the 
geometrical symmetry of C50CI10 where it has the D5/1 sym- 
metry in xy plane with diameter (10.3 4) almost twice larger 
than that along z axis (5.4 A), we find the crystal lattices such 
as tetragonal (tt), simple hexagonal (sh), fee and hep struc- 
tures could be the most possible candidates. The former three 
are shown in Fig. 2, where the coordination systems of the lat- 
tices are coincident with those of the C50CI10 molecule. Then, 
the total energy is calculated accordingly with respect to dif- 
ferent lattice parameters for each presumed structure L2QJ. The 
total energy is scanned in every OA A for lattice constants a 
and c. Owing to the limits of the method and the capability of 
our computer system, a properly fixed orientation of C50CI10 
molecules, namely, the xy plane of C50CI10 molecule is ar- 
ranged to be vertical to the z axis in tt, sh, and hep structures, 
and to be along the (11 1) plane in the fee structure, is assumed 
in calculations, which makes the system bear a relatively high 
symmetry fsT]. 

Table I gives the calculated minimal relative total energies 
per C50CI10 molecule and the corresponding lattice constants 
for the four presumed solid phases, where the data with pre- 
cision of 0.01 A is obtained by a bicubic interpolation from 
the calculated data with precision of 0. 1 A. For both tt and sh 
structures, the total energy shows a series of minimums (la- 
beled by A-F) with almost the same values, suggesting that 



these two structures are labile. For the hep and fee structures, 
only one minimal total energy is obtained, showing that the 
hep and fee C50CI10 could be formed, though the hep phase 
looks more stable than the fee phase in energy per molecule, 
as the minimal energy of the hep structure is about 0.6 eV per 
C50CI10 molecule lower than those of tt, sh and fee structures. 

To show it clearly, the contour map of the total energy for 
the hep C50CI10 in the plane of lattice constants a and c is 
depicted in Fig. 3(a), where only one minimum of the total 
energy is observed at a = oq =13.57 A, and c — cq =11.39 
A, as marked by the letter A. Figs. 3(b) and (c) present the 
total energy per C50CI10 molecule of the hep C50CI10 crystal, 
which is measured from the total energy of isolated C50CI10 
clusters 1 32], as a function of lattice constants a and c, re- 
spectively. The results reveal that the hep structure is stable. 
Fig. 3(d) gives the total energy per C50CI10 cluster in the fee 
C50CI10 crystal as a function of lattice constant a, where only 
one minimum of the total energy is seen at a = aq =19.20 
A, indicating that the fee C50CI10 crystal could also be stable, 
though the minimal total energy of the fee phase is larger than 
that of the hep phase. 

In order to confirm the relative stabilities of the obtained 
hep and fee lattices, we have optimized the atomic positions 
with fixed lattice parameters by 40 steps Broyden-Fletcher- 
Goldfarb-Shanno minimization. The results show that no no- 
table atomic displacements are observed, and the total-energy 
changes are less than 0.0003% (0.04eV), indicating that the 
C50CI10 molecules keep their shapes and sizes, and the lattice 
structures which we have obtained are stable. 

The calculated cohesive energy per C50CI10 molecule is 0.9 
eV for the hep solid, and 0.3 eV for the fee solid, which is 
much smaller than the typical C-C bond energy (more than 3 
eV). The cohesive energies of both solid phases are less than 
that offcc Ceo solid (1.6 eV^) 1 20], implying that C50CI10 solid 
is less stable than Ceo solid. The contour map of the valence- 
electron density for hep C50CI10 crystal is given in Fig. 4. 



TABLE I: The calculated minimal relative total energies per C50CI10 
molecule and the corresponding lattice constants for the presumed 
four solid structures. 







Lattice constants (A) 
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Energy (eV) 
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13.99 


9.22 


0.77 




B 


13.61 


9.72 


0.77 


tt 


C 


14.02 


9.80 


0.76 




D 


13.72 


10.31 


0.75 
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13.99 


10.82 


0.73 




F 


13.65 


11.41 


0.74 




A 


13.37 


9.75 


0.55 




B 


13.38 


10.42 


0.58 


sh 


C 


13.33 


10.79 


0.57 
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13.39 


11.36 


0.54 
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13.38 


12.37 


0.54 




F 


13.38 


12.99 


0.55 


hep 


A 


13.57 


11.39 


0.00 


fee 


A 


19.20 




0.60 



It is seen that the valence-electron density between C50CI10 
molecules is much lower than the density in a single C50CI10 
molecule. The similar situation occurs for the fee phase of 
C50CI10. 

The energy bands, DOS G{e), as well as the integrated 
DOS I{E) = J^^ G{e)de for solid C50CI10 are calculated 
for the hep and fee phases with the optimized lattice constants 
flo and cq. Fig. 5 presents the band structures of the hep and 
fee C50CI10 crystals, respectively, where only several energy 
bands around the energy gap are shown. For the hep crystal, 
the energy difference between the valence-band top (VBT) at 
the r point and the conduction-band bottom (CBB) at the K 
point is 1 .59 eV, and for the fee crystal, the energy gap be- 
tween the VBT at the F point and the CBB at the L point is 
about 1.79 eV, showing that both hep and fee C50CI10 solids 
may be indirect-gap band insulators. The present situation is 
in contrast to the fee Cgo solid which is a direct-gap insula- 
tor with an energy gap of 1.5 ey|20|. We have also calcu- 
lated the energy levels of C50CI10 molecule, and found that 
the energy gap between the highest-occupied molecular or- 
bital (HOMO) and the lowest-unoccupied molecular orbital 
(LUMO) is around 1 .9 eVVwhich is in well agreement with a 
previous calculated resultf 1 Ij, and is comparable with that of 
Cgo cluster where the energy gap between the HOMO state, 
hu, and the LUMO state, is about 1.9 6^120]. Obviously, 
the energy gap remains finite when the C50CI10 molecules 
are condensed to form a solid. Like other typical molecu- 
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FIG. 3: (Color online) (a) The contour map of total energy for hep 
C50CI10 in the plane of lattice constants a and c, where A indi- 
cates the minimum of the total energy. The total energy per C50CI10 
molecule, measured from the total energy of the isolated C50CI10 
molecule, as a function of the lattice constant a (b) and c (c) for hep 
phase, where the optimized lattice constants ao = 13.57 A, co ~ 
1 1 .39 A, and (d) for fee phase with the optimized lattice constant 
ao = 19.20 A. 




FIG. 4: (Color online) Contour map of the valence-electron density 
of hep C50CI10 crystal in the (001) plane. The larger lobe-like pat- 
terns stand for those of CI atoms, while the inner smaller dumbbell- 
like patterns stand for those of C atoms, where the highest-density 
contour (red, around C and CI atoms) represents the equal-density 
line of 1.383 M^, and the lowest-density contour (blue, dark grey) 
represents the equal-density line of 0.054 /A^. It shows that the 
interactions between the CI atoms in the two neighboring C50CI10 
molecules are quite weak. 
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FIG. 5: (Color online) Band structures of hep (a) and fee (b) C50CI10 
around the energy gap, where the Fermi level is set to zero. For 
the hep phase, the valence-band top (VBT) is at F point, and the 
conduction-band bottom (CBB) is at K point; for the fee phase, the 
VBT is at the F point, and the CBB is at the L point, showing the 
solid C50CI10 is an indirect band insulator with a gap of 1.59 eV 
(hep) or 1.79 eV (fee). 



lar svstems ll33ll . the energy bands of these two C50CI10 crys- 
tals become narrow, as they are formed through the overlap of 
C50CI10 molecular energy levels. 

Fig. 6 presents the DOS and the integrated DOS of the 
hep C50CI10 as an example, which are obtained by means 
of a smearing technique. An energy gap in the vicinity of 
the Fermi level is clearly seen. From the positions of peaks 
in the DOS around the energy gap, the information of the 
HOMO-derived and LUMO-derived bands as well as next 
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higher bands can be extracted. Evidently, the DOS of elec- 
trons at the LUMO-derived states appears to be smaller than 
that at the HOMO-derived states. With the increase of en- 
ergy, the integrated DOS increases nonmonotonically below 
and above the energy gap, where the flat part in the vicinity 
of the Fermi level measures the width of the energy gap. The 
presented results in Fig. 6 are in agreement with the calcu- 
lated electronic energy levels shown in Fig. 5. It should be re- 
marked that the obtained DOS as a function of energy shows a 
very similar behavior for the hep and fee phases, though there 
are slightly quantitative changes. The reason is that owing 
to the quite weak interactions between C50CI10 molecules, as 
manifested in Fig. 4, the DOS of electrons for solids would 
have profiles similar to that for clusters, yielding that the DOS 
is similar for both hep and fee C50CI10. 

Let us now address another interesting issue if C50CI10 
molecules can condense to form a monolayer structure. Re- 
call that Ceo clusters can form a hexagonal monolayer struc- 
ture which has gained much attention recently! 34 1. From Ta- 
ble I, one may note that for the sh lattice, the minimal total 
energies have slight changes with the lattice constant a, but 
have remarkable changes with constant c. It hints that a sta- 
ble layered structure with a single minimal total energy could 
probably be formed. To probe this possibility, we have taken 
the lattice constant c of the sh supercell to 40 A to decline 
the interactions between the C50CI10 layers but keeping the 
molecular orientation the same as in the sh lattice. It turns out, 
after some efforts, that a single minimal total energy with the 
optimized lattice constant aq = 13.52 A is indeed obtained. 
The cohesive energy per C50CI10 molecule for this monolayer 
is found to be 0.3 eV. The energy bands and DOS show that 
the C50CI10 hexagonal monolayer is a band insulator with an 
energy gap of 1.85 eV. 

It would be emphasized that as manifested in Fig. 4, 



the molecule-molecule interactions in the solid C50CI10 ap- 
pear to be weak. Whether the interactions between C50CI10 
molecules in the solid phase are of the van der Waals force 
or not, is at present unclear within our calculations, because 
it is known that the LDA cannot correctly describe the van 
der Waals force as it assumes fundamentally that the electron 
density is flat like in metals (see, e.g. |35]). Nevertheless, our 
calculations indeed suggest that the molecular solid C50CI10 
with either the hep, fee or hexagonal monolayer phase can be 
somehow formed through weak interactions. 

To summarize, in terms of first-principles calculation 
within the density functional theory, we have predicted a 
molecular solid C50CI10 with either a hep, fee, or a hexag- 
onal monolayer structure. These possibly stable structures are 
determined from the total-energy calculations, where the hep 
phase is found more stable than the fee phase and the hexag- 
onal monolayer in energy per molecule. The energy bands 
and DOS for the hep and fee C50CI10 as well as the hexago- 
nal monolayer are discussed. Our results show that C50CI10 
solid can be a band insulator with an indirect energy gap of 
1.59 eV for the hep structure and 1.79 eV for the fee struc- 
ture, or with a gap of 1.85 eV for the hexagonal monolayer. 
Since the LDA calculation usually underestimates the energy 
gap iiH Isiil . our present simulations inevitably call for fur- 
ther theoretical studies ll37ll or experiments. We expect that 
the C50CI10 solids can be successfully synthesized in near fu- 
ture. 

The aufliors are grateful to B. Gu, H. F. Mu, X. B. Wang, 
and Z. C. Wang for helpful discussions. This work is sup- 
ported in part by the National Science Foundation of China 
(Grant Nos. 90403036, 20490210), and by the MOST of 
China (Grant No. 2006CB601 102). 
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FIG. 6: (Color online) Density of states (DOS) and the integrated 
DOS (dashed line) of electrons as a function of energy for the hep 
C50CI10. The energy zero-point is taken at the Fermi level, that is 
indicated by the vertical dash-dot line. The flat part of the integrated 
DOS in the vicinity of the Fermi level gives the width of energy gap. 
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